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A counting formula for computing the number of (Majorana) zero modes bound to topological 
point defects is evaluated in a gradient expansion for systems with charge-conjugation symmetry. 
This semi-classical counting of zero modes is applied to some examples that include graphene and a 
chiral p-wave superconductor in two-dimensional space. In all cases, we explicitly relate the counting 
of zero modes to Chern numbers. 



I. INTRODUCTION 

The counting of zero modes, cigenstatcs annihilated 
by a single-particle Hamiltonian ri, has a long history 
in physics. Charge-conjugation symmetry, the existence 
of a norm-preserving linear (antilinear) transformation 
C that anticommutes with %, protects the parity of the 
number of zero modes. When the parity is odd, at least 
one zero mode must be robust to any perturbation that 
preserves the charge-conjugation symmetry. This pa- 
per aims at calculating the parity of zero modes of a 
single-particle Hamiltonian H(p, <p{xj) when (1) it obeys 
charge-conjugation symmetry, (2) it describes fermionic 
quasiparticles, and (3) it depends on a position depen- 
dent vector- valued order parameter cp(x). These three 
assumptions are often met in mean-field treatments of 
electrons interacting with each other or with collective 
excitations such as phonons or magnons in condensed 
matter physics. 

Zero mode solutions can be found by direct means, in 
practice solving a differential equation. This requires a 
non-universal definition of the model since both micro- 
scopic and macroscopic data must be supplied, say the 
boundary conditions to be obeyed at the origin and at 
infinity in space. 

Is there an alternative approach to calculating the par- 
ity in the number of zero modes that is more universal? 
The celebrated index theorem for elliptic differential op- 
erators gives a positive answer to this question for those 
problems in physics for which this theorem applies (Dirac 
Hamiltonians for example) .El The index theorem achieves 
this by relating some (not all!) zero modes to a topolog- 
ical number (a global property of the Hamiltonian that 
can only take integer values). However, the index theo- 
rem cannot be applied to most Hamiltonians of relevance 
to condensed matter physics. 

In this paper, we start from an exact integral repre- 
sentation of the total number of unoccupied zero modes 
(up to exponential accuracy) N for fermionic single- 
particle Hamiltonians with charge-conjugation symmetry 
in terms of a conserved quasiparticle charge Q, 

N = -2Q. (1.1) 
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FIG. 1: The charge-conjugation-symmetry-breaking param- 
eter <f> moves in energy a mid-gap state upward or downward 
depending on its sign. The continuum parts of the energy- 
eigenvalue spectrum are denoted by the shaded boxes. The 
thin horizontal line denotes the band center about which the 
spectrum is symmetric when <fi = 0. 

This counting formula appears implicitly in Rcf. ^ and 
explicitly in Ref. ||[ both in the context of polyacetylene.cl 
For polyacetylene, it relates the number of unoccupied 
zero modes to the conserved electric charge Q dw induced 
by domain walls in the spontaneous bond ordering trig- 
gered by the coupling of electrons to phonons. Remark- 
ably, the electric charge Q dw = ±1/2 induced by a single 
domain wall is fractional and counts a single zero mode, 
with the sign ambiguity jresolved by whether the midgap 
state is filled or empty.&Lj Alternatively, this sign ambigu- 
ity can be removed by the application of a small charge- 
conjugation-symmetry-breaking perturbation that shifts 
the energy of the zero mode up or down (see Fig. [l]). 
Hence, the counting formula ([O]) becomes 

N = 2Q mod 2 (1.2) 

if no prescription is given as to whether the filled Fermi 
sea includes or not a zero mode. The same assignment 
of quantum numbers also relates a single zero mode and 
the conserved electric charge Q Kcku i c induced by a vortex 
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with unit j/octicity in the Kckulc dimcrization pattern of 
grapheneJxEj 

For polyacetylene and graphene, the charge- 
conjugation symmetry is approximate, for it origi- 
nates from a sublatticc symmetry that is broken as 
soon as next-nearest-neighbor hopping is included 
in the tight-binding model. To the extend that the 
Bardccn-Cooper-Schrieffcr (BCS) mean-held approx- 
imation to superconductivity is empirically observed 
to be excellent, single-particle Bogoliubov-de-Gennes 
(BdG) Hamiltonians realize a much more robust charge- 
conjugation symmetry.E3 For BdG Hamiltonians, zero 
modes are associated to Majorana fermions. Majorana 
fermions do not carry a well-defined electric charge since 
global electro-magnetic gauge invariance is broken in 
any mean-held tr eatm ent of superconductivity. The 
counting formula (1.1) nevertheless applies to any BdG 
Hamiltonian H BdG with the important caveat that the 
conserved quasiparticle charge <3 BdG is unrelated to the 
electric charge. Rather, the conservation of QbcIG cn ~ 
codes, for any local BdG Hamiltonian, a local continuity 
equation obeyed by the Bogoliubov quasiparticles that 
is responsible for the conservation of the thermalpflpw in 
the mean-field treatment of superconductivity£irli3 



the space coordinate 



In this paper, we represent the counting formula (1.1) 
in terms of single-particle Green functions. The advan- 
tage of this choice is that it easily lends itself to a pertur- 
bative (gradient) expansion of the conserved quasiparti- 
cle charge. Q that forgoes the non- universal short-distance 
dataJiZHiia To leading order, this expansion is akin to an 
adiabatic approximation. We thus propose the adiabatic 
approximation to the conserved quasiparticle charge Q 
as an efficient mean to compute the parity in the num- 
ber N of unoccupied zero modes of charge-conjugation- 
symmetric single-particle Hamiltonians. 

We then apply this formula to charge-conjugation- 
symmetric Dirac insulators with time-reversal symmetry 
in arbitrary dimensions and chiral p-wave BdG super- 
conductor in two-dimensional space that all support a 
point defect. One of the main results of this paper is an 
integral representation for the number (zero or one) of 
unoccupied zero mode induced by a unit vortex in a two- 
dimensional chiral p-wave BdG superconductor. A by- 
product of this integral representation is that it is closely 
related to the second Chern number. We also relate the 
d-th Chern number to the number of unoccupied zero 
modes for Dirac fermions in d-dimensional space that in- 
teract with a d-tuplet of Higgs field supporting point de- 
fects. 

The gradient expansion is pre sented in Sec. [n] and ap 



plied to point defects in Sec. III. We conclude in Sec. IV 



II. 



GRADIENT EXPANSION OF THE 
COUNTING FORMULA 



Let the charge-conjugation-symmctric single-particle 
Hamiltonian H(p 7 <p(x)) be such that its dependence on 



x e 



(2.1) 



is implicit through that of a static vector-valued order 
parameter 

<p(x) E R D 7 (2.2) 
while its dependence on the canonical momentum 



P 



id EE -i — 

dx 



(2.3) 



is explicit. (We have chosen natural units, i.e., H = 1.) 

The Hermiticity of the single-particle Hamiltonian 
H(p,<p(x)) implies the existence of a conserved quasi- 
particle charge, which we construct explicitly with the 
help of the second quantized HamiltonianEJ 



Hee j &(x)H(p,ip(x))%(x). 



(2.4a) 



The creation operators ty'(x) and the annihilation oper- 
ators ^(x) obey here the fermion algebra 



[%(x) t %(x')}=0, r,r' = !,••• ,R, 



(2.4b) 



in some i?-dimensional representation of the charge- 
conjugation-symmetric single-particle Hamiltonian 



H(p,<p(x)) = - C- 1 H(p, V (a;))C. 



(2.4c) 



The norm-preserving operation for charge conjugation 
in the i?-dimensional single-particle representation is de- 
noted by C. The short-hand notation J stands for the 

X 

space integration J d d x. The definition of the conserved 
quasiparticle charge requires filling of the Fermi sea for 
the quasiparticles created by the field ^l(x). 

We begin with an expression for a quasiparticle charge 
density in the vicinity of a point x. To do so, we shall 
assume that the static vector- valued order parameter can 
be decomposed additively, i.e.JEl 



ip{x) = ip + Sip{x), 



(2.5) 



in such a way that: (i) the single-particle Hamiltonian 

H q (p)eeH{p, Vq ) (2.6) 

is translation invariant, (ii) the single-particle eigenvalue 
spectrum of H (p) is fully gaped with the gap 2A > 0, 
(hi) the changes between 6<p(x) and 5(p(y) when x and y 
are a distance of order 1/A apart are small. Condition 
(iii) implies that the gradient of Sip(x) can be viewed 
as a smooth and small perturbation that can be treated 
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perturbatively in the gradient expansion that will follow 
shortly. 

Condition (i) implies that 



(2.7a) 



with the symmetric Fourier conventions 

*t(sc) := (2n) +d/2 J e~ ip:c &(p), 
p 

§t(p) := (2n)- d / 2 J e +i r m $\x), 



(2.7b) 



for the annihilation (creation) operators and the asym- 
metric convention 



K(x) := J e +ip x K,(p), 
p 

K(p) := / e~ ipx K(x), 



(2.7c) 



for any kernel /C (such as the Hamiltonian %). Here, J is 

p 

a short-hand notation for the momentum-space integra- 
tion / d d p/(2n) d . 

Condition (ii) implies the existence of the characteris- 
tic length scale 1/A . 

Define the quasiparticle charge density 

/j 
— {x\tx R [g{u)-g o {u)]\x). (2.8a) 



Here, tr R denotes the trace over the i?-dimcnsional de- 
grees of freedom and we have introduced the Euclidean 
single-particle Green functions 



i 



iui-W 



1 



iu> — H a 



(2.8b) 



The quasiparticle charge density p 7 (x) depends on the 
contour of integration 7. The latter is chosen as in Ap- 
pendix |A| so that the integration over the cj-complex 
plane picks up only the first- order poles from the non- 
vanishing an d ne gative energy eigenvalues of % and of 
n [see Eq. (|A15|) 1. 

It is s hown in Appendix |a] [see Eq. ( A16c )] that 
Eq. ( |2.4cD and condition (ii) imply that the total num- 
ber N of unoccupied zero modes of % is related to the 
conserved quasiparticle charge Q 1 by 



iV 
~2 



- = - I d d xp y (x) = -Q r 



(2.9) 



According to the counting formula (2.9), computing N 
reduces to computing the conserved quasiparticle charge 
Q 1 induced by the s moo th variation of Stp(x) through 
space defined in Eq. (2.1). The Hcrmiticity of the single- 
particle Hamiltonian (2.4c) guarantees the existence of 
this conserved quasiparticle charge. A local law for the 
conservation of quasiparticle current follows when the 
single-particle Hamiltonian is local. For single-particle 
Hamiltonians with the global U(l) symmetry delivering 
the conservation of the total electric charge, the con- 
served quasiparticle charge Q 7 is nothing but the elec- 
tric charge in units in which the electron charge is unity. 
For single-particle Hamiltonians describing Bogoliubov- 
de-Gennes (BdG) quasiparticles, this global U(l) symme- 
try is spontaneously broken. The conserved quasiparticle 
charge Q 1 is then related to the conserved thermal cur- 
rent of BdG quasiparticles. The counting formula 
appears implicitly in Ref. ^| and explicitly in Ref. |^.C 

There is an alternative to specifying 7 in the count- 
ing formula (2.9). We can regulate the first-order pole 
of the Green function G(uj) at us — by adding a per- 
turbation that moves all zero modes to strictly positive 
energies. This perturbation must be small if all these 
positive energy increments are to remain much smaller 
than the threshold A to the continuum. We can then 
safely replac e th e contour of integration 7 in the count- 
ing formula ( [2.9[ ) by R since the subtraction of Q from 
Q insures the convergence of the us integration for large 
ui. 

For example, we imagine that it is possible to aug- 
ment the vector- valued order parameter (2.5) by the 
conjugation-symmetry-breaking real-valued field <fi with- 
out loosing conditions (i)-(iii). More precisely, we define 
the (D + l)-tuplet 



00*0 







( fl( x )\ 


(j} D { X ) 




ip D (x) 






{ * J 



(2.10a) 



and we assume that Eq. ( 2.4cj ) becomes 

n(p,<p(x),<i>) = -c- 1 H{p,t P (x),-4>)c. 



It then follows that 



Q{4>) 



/ d " xp(x " 



-Q{- 



(2.10b) 



(2.11) 



where the quasi part icle charge density p(x, <p) is ob- 
tained f rom E q. (2.8) with R substituting for 7, Hamil- 
tonian (2.10b) substituting for H(p, ip(xf\, and T~L = 
H{p,(Pq,4>) substituting for Ti a = H(p,<4>o)- The smok- 
ing gun for the unoccupied zero modes is now a discon- 
tinuity at the conjugation-symmetric point = o f the 
odd function Q(<j)) of </>, i.e., the counting formula ( |2.9| ) 
has become 



N 

— = — lim d x p(x, ( 
2 *-*>.' 



- lim Q{<t>) 

0-s-O 



(2.12) 
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where the sign of <f> is to be chosen so as to move the 
N zero modes along the energy axis to positive energies. 
Furthermore, we can relax the condition that the sym- 
metry breaking <j) is constant everywhere in space (x) 
provided that the condition 



fa 2 (x) 



A 2 



(2.13) 



holds everywhere in R d . Condition (iii) then means that 
5<p{x) varies slowly on the characteristic length scale 
1/A . 

Rather than computing p{x) exactly, say with the help 
of numerical tools, we are after the leading contribution 
to the gradient expansion of the quasiparticle charge den- 
sity p(x), which we shall denote as /f adia (a;) where the 
subscript "adia" refers to the adiabatic approximation 
contained in condition (iii). 



The order parameter (2.10a) enters linearly in all the 
single-particle Hamiltonians that we shall consider ex- 
plicitly in this paper. Hence, there follows the additive 
law 

n(p, 0(a)) =U(p, fa ) + V(p, Sfa(x)) 
= H + V(p,Sfa(x)) 

upon insertion of 

0(x) = fa + Sfa(x). (2.15) 
Upon second-quantization, this implies that 

H=J J&(p) [H (p)6(p-q)+V(p;q)]$(q) 
p i 



(2.16a) 



holds with 



V(P;9):= (f) (^)-^- 9) " (2 ' 16b) 

The subscript means setting Sfa to zero so that the 
gradient 



dcj)) V dep , 



(2.16c) 



depends only on the single-particle canonical momentum 
operator p = —id. We have also adopted the convention 
that matrix elements of p are to be symmetrized, i.e., 

v*>w--i(r (||)-(£),)(.> 9-m 

for any differentiable functions / and g. 

We now expand the quasiparticle charge density up to 
the first non-trivial order in an expansion in powers of V 
with the help of the geometrical series 



0(w)-0 o H«£(0oMvra o M 



(2.17) 



This gives the expansion 

oo 

p(x) *iJ2p n (x) 
n=l 

with 

Pn( X ) 



(2.18a) 



= i(<?iH hq„)-a; 



" P ?i q„ 

D+l 

X 

a, , ■ ■ ■ ,a„ — 1 



/ a„,-,a 1 ( w ,P:<Zl,' 

x S( P an {q n ) ■ ■ ■ H 3i (qi) 

where / = / ^ and the integrand 

4„,.., ai ( w ,p>9i>'" ,q n ) 



• In) 



(2.18b) 



:=tr. 



£o( w >P + <Zl+<?2 + ••• + <?«) 

an x 



\ " / u 

x Q ( w ,P + <?i +Q2 

( on 



On 

2 



90a 



■+Qn-l) 

u},p + q 1 +q 2 H hq„_ 2 



on 



ui,p - 



9i 



G {u,p) 



(2.18c) 



Finally, we expand the integrand (2.18c) in powers of 
the coordinates ■■• , q ni of the momenta q 1 , ■ • • , 
q n , which are then to be integrated over. Condition (iii) 
suggests that we keep only first-order derivatives in the 
slow ly v arying fluctuations 8<p(x) of the order parame- 
ter (2.5). This approximation yields the leading contri- 
bution p adla (x) in the gradie nt expansion of the quasi- 
particle charge density ( 2.18a ) and becomes exact in the 
limit when the ratio between the characteristic length 
scale 1 / A and the characteristic length scale over which 
the change in 6<p(x) becomes significant vanishes. 

However, to each order n in this expansion there are 
terms for which not all Sfa(x) are differentiated. These 
terms do not have to be small. Hence, they should be 
treated non-perturbatively. This is achieved by replacing 
the Euclidean single-particle Green function 



G (u,p) 



1 



Hip, fa 



with the Euclidean semi-classical Green function 

1 



G^ c (u,P,x) := : 



n=l 



ho - H(p,fa(x)) 



(2.19) 



(2.20) 
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to any given order n in the expansion (2.18). This result 
is rooted in the fact that the choice of 4> in the additive 
decomposition ( 2.15 ]) is arbitrary from the point of view 



of the expansion ( 2.18 ). Such arbitrariness should not 
appear in the final result, i.e., the final result should only 
depend on 4>(x) orJL<f) a (x) = ^r(x) with a = 1, • • ■ ,d+ 
1 and i = 1, • • • , dS3 

Thus, by collecting the appropriate derivatives of the 
spatially varying order parameter 5<fi(x), we arrive at 
expressions for the induced quasiparticle charge density. 
We now analyze this expression according to which term 
is the first non-vanishing contribution to the expansion. 

When the adiabatic approximation for the quasiparti- 
cle charge density is the non- vanishing n = 1 term in the 
gradient expansion, it is given by 

PadiaN=\a 1 (-i)(^ 1 ai ) (*). (2-21&) 

Here, the summation convention is assumed over re- 
peated indices and the coefficients are given by 



1 



:tr, 



U) p 



dg- 1 dg 



g\ (w,p). (2.21b) 

/ 



The subsc ript refers to the semi-classical Green func- 
tion ( 2.20 ) . This case is the relevant one for the study of 
point defects in one-dimensional space. 

When the first non-vanishing contribution to the adi- 
abatic expansion occurs at n = 2, then the quasiparticle 
charge density contains two gradients and is given by 



PadiaW 



iiai (-i) 2 (^>a 2 ) (9iA)0»0- ( 2 - 22a ) 



By assumption, X i a defined by Eq. (2.21b) vanishes, but 



X 




dg dg-'dg- 1 dg 

2- — — y- 



dp t d(j), d<j>, dpi 



-Q- — 7T— Q (w,p) 



(2.22b) 



does not. Again, the summation convention is assumed 
over repeated indices and the subs cript refers to the 
semi-classical Green function ( 2.20 ). We have chosen to 
represent Eq. (2.22) so as to make the reality of p^J^x) 
manifest. This case is the relevant one for the study of 
point defects in two-dimensional space. 

Observe that, whenever n > 1, we must allow for the 
possibility that 



d 2 n 

9Pid(j) a 



d'g- 1 

d Pi d4> a 



(2.23) 



is non- vanishing for some i = 1, ••■ ,d and some a = 
1, •••,£) + 1. These terms occur when dealing with a 
p-wa vc su perconductor in d = 2 dimensions as we do in 
Sec. IIIB ; however, we find by explicit calculation that 
these terms vanish upon integration over u and p. 

When the order parameter is independent of momen- 
tum, 



d 2 n 

dp t d4> a 



(2.24a) 



for any i — 1 , • ■ • , d and any a = 1, • • • , D + 1. For 
example this is the case with D = d for Dirac fcrmions in 
d-dimensional space interacting with (d + 1) jccal- valued 
Higgs fields, in which case it is the coefficients 



tr, 



, d Pt 30 



dp. 



duj 



'(2.24b) 



that controls the adiabatic approximation to n-th order 
through 

PadiaW = (-i)% n3n ... liai (S in ^ a „) ' ' ' (\4> a J ( x ) 

(2.24c) 



as we shall show in Sec. IIIC. Th e sub script refers to 



the semi-classical Green function (2.2C). 

Finally, with the expression for the local quasiparti- 
cle charge density p a( ji a ( a; ) within the adiabatic approxi- 
mation in hand, we can compute the total quasiparticle 
charge Q adia - Naturally, one goes from the local den- 
sity to the total charge by integrating the former over all 
space. We conclude that 



(2.25a) 



Qadia — / ^ x PadiaC 31 ) 

= J d d x(-i) d l inan ... ii3i (x) 

x (d i J an )---(d ii <f> ai )(x), 

where it is the Euclidean semi-classical Green function 

1 



S s _ c (w,p,a;) 



that enters the kernel X. 



-H(p,<f>(x)) 



(2.25b) 



III. ZERO MODES INDUCED BY POINT 
DEFECTS 

We are going to apply the adiabatic expansion from 
Sec. H to the case of point defects supported by the static 
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order parameter 



Generic single-particle H when d — 1 



<f>(x) = 







( 


4> D (x) 




ip D (x) 


\<t>D+l{v)j 




{ 4>(x) J 



pd+1 



(3.1a) 



and 



4>\x) 



Vo 



when space is d-dimensional, i.e., when 
x e R d . 



(3.1b) 



(3.1c) 



The component (f> D+1 (x) = <fr(x) breaks locally the 
charge conjugation symmetry. This component deter- 
mines if a zero mode bound to a defect at x is occu- 
pied or unoccupied. All remaining components of the 
order parameter (3.1a) are compatible with the charge- 
conjugation symmetry of the single-particle Hamiltonian. 
The condition (3.1t) suggests that the homotopy group 



Tl d (S d ) = Z 



(3.2) 



of smooth maps from the compactification of M. d into the 
d-sphere S d to the d-sphere S d in ostLer-parameter space 
K d+1 might play an important role.E3 

We begin in one dimensional space (x) with a generic 
single-particle Hamiltonian. We show that the conserved 
quasiparticlc charge Q adia that enters the counting for- 
mula (^^) is the first Chern number if we relax the 
condition that the charge-symmetry-breaking component 
4> of the order paramet er is infinitesimally small and 
if we compactify space (3.1c) and the order-parameter 



space (3.1a 



We continue with the chiral p-wave superconductor in 
two-dimensions when the superconducting order param- 
eter supports a vortex. We show that the number of 
unoccupied zero modes bound to a vortex with unit vor- 
ticity computed from the adiabatic approximation agrees 
with the direct construction of zero modes once all mi- 
croscopic data have been supplied. Moreover, we show 
that the conserved quasiparticlc charge Q adia that enters 



the counting formula (2.9) is also related to the s econd 
Chern number after compac tifica tion of both space (3.1c) 
and order-parameter space ( 3.1a ). This is apriori surpris- 
ing since the superconducting order parameter couples 
to the momentum contrary to the simpler case of Dirac 
fermions. 

We also study Di rac si ngle-particle Hamiltonians in d- 
dimensional space ( 3.1c ) that are represented by Dirac 
matrices of dimension R = 2 d . We show how the con- 
served qu asip article charge Q adia that enters the counting 
formula ( |2.9| ) is related to the d-th Chern number when 
the Dirac spinor couples to a (d+ l)-tuplct of Higgs field. 

We close by discussing how to interpret the adiabatic 
approximation. 



We compactify space and the order-parameter space, 
x g S\ (f)(9) e S 1 C M 2 , (3.3) 
respectively, such that 



<£(x) 



=m ( sin0(x) 
'x)l lcos6>(x) 



J 2 



(3.4) 



and denote with Q = (icj — TL)^ 1 the single-particle Green 
function in Euclidean space for any suitable RxR matrix- 
valued single-particle Hamiltonian Hvp, 0(x)J. Suitabil- 
ity means that the dependence on the momentum oper- 
ator p = —id x in H(p, </>(x)) is only restricted by locality 
while that on the two-component order param eter c b(x) 
is linear. Furthermore, H(p,<f>(x)) obeys Eq. ( |2.10| ) for 
D = 1 and supports the spectral gap 2m for a uniform 
4> . Under these conditions, for any p the semi-classical 
RxR matrix H (p, 4>(x)) is traceless and its square is pro- 
portional to the unit matrix with the smallest nonvanish- 
ing eigenvalue no smaller than to 2 . We have introduced 
the spherical coordinate 8 of the one-sphere S 1 C R 2 in 
order-parameter space. 



charge Q adia becomes 



Qadia — 



According to Eq. (2.21), the conserved quasiparticle 



-tr R (gd p g- l Qd s Q- l g 



w pes 1 xes 1 



(3.5) 



GdeG-^g-^y^dj. 



The subscript refers to the semi-classical Green func- 
tion (2.20). We use the identity 



(3.6) 



d^ 1 = i, 

introduce the family of indices labeled by v, 

^1)^2)^3 = MA (3-7) 
and the Euclidean 3-momentum 

(u,p,6). (3.8) 



With the help of t he manipulations made bet ween 
Eqs. (|C3j) and (|C45|) , it is possible to re- write Eq. (|J) 



as 



2tt 



2 71 



^adia 



24tt 2 



doj / dp / dO e 



(3.9) 



xtrjGdvG-'GdvG-'GdvG- 1 



Equation ( |3.9| ) is the first Chern number^ Thus, the con- 
served quasiparticle charge Q aA - obey ing Eq. (3.9) t akes 
integer values. We defer to Sec. HID and Appendix CI 
for a more detailed discussion of the connection between 
the first Chern number and the number of unoccupied 
zero modes. 
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B. Chiral p-wave superconductor when d = 2 

1. Definition 

One of the main results of thi s paper consists in apply- 
ing the counting formula (2.12) to the two-dimensional 
chiral p-wave superconductor with the single-particle 
BdG Hamiltonian 



H 



BdG 



\{Pi 

-ip 2 ,A*(x)} 



ip 2 ,A(aj)} 
-e(p) 



BdGf 



P_ 

2m 



e{p) 

(3.10) 

for the case when the superconducting order parameter 
supports the vorticity n v = ±1 at the origin of two- 
dimensional Euclidean space, i.e., when 



A(x) :=A (r)e"^ 
where r and 9 are the polar coordinates of x 6 T 



(3.11) 



The applicability of the counting formula (2.12) follows 
from the antiunitary conjugation symmetry 
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The Pauli matrices p ly p 2 , and p 3 encode the particle and 
hole block structure of the single-particle BdG Hamilto- 
nian (3.10). Complex conjugation is denoted by K. 

For a unit vorticity, Read and Green have shown in 
Rcf. ^6| the existence of a single zero mode bound to the 
vortex by solving the eigenvalue problem at zero energy. 
A zero mode in a single-particle BdG Hamiltonian real- 
izes a Majorana fermion, for it cannot be distinguished 
from its particle-hole conjugate. For an arbitrary vor- 
ticity n v G Z , it is shown in Refs. |27] and |2^ by solving 
the eigenvalue problem at zero energy that the number of 
Majorana fermions is one if n v is odd and zero otherwise. 

We do not expe ct th e adiabatic approximation to the 
counting formula ( 2.12] ) to capture this subtle parity ef- 
fect since it is only sensitive to the net vorticity n v 
trapped in region of space much larger than the char- 
acteristic length scale I 3> 1/A over which |<5tp| changes 
by the amount of order A . The adiabatic approxima- 
tion fails to distinguish the case of a single vortex with 
vorticity n v and n v vortices with unit vorticity that are 
separated by a distance of order i. This parity effect is 
a non-pcrturbativc effect from the point of view of the 
gradient expansion that we now present. 



2. Counting zero modes 

To count the unoccupied zero modes induced by a vor- 
tex in the superconducting order parameter with the gra- 



dient expansion of Sec. |l], we need to simultaneously 
move in energy the zero mode and properly regulate the 
singularity at the core of the vortex. Achieving both 
goals is impossible with the 2x2 BdG Hamiltonian for 
the chiral p-wave single-particle Hamiltonian. On the 
other hand, both conditions are met after doubling the 
BdG single-particle Hamiltonian so as to introduce an 
auxiliary chiral-symmctry-breaking perturbation in two 
ways (denoted by the subscripts =f), 

^<**<»>>4i«. {££<»>> -«H, } 

(3.13a) 

The short-hand notations := p x ± ip y was introduced 
and p is the 2x2 unit matrix. The triplet 
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is real-valued and is made of the static superconduct- 
ing order parameter A(x) = A 1 (a;) + iA 2 (a;) and of the 
static auxiliary charge-conjugation-symmetry-breaking 
field (j) 3 {x). 

The spectru m of t he pair of auxiliary single-particle 
Hamiltonia ns ( 3T3a) given a uniform order parameter 
O in Eq. ( ft.l3b| ) is 



(3.14) 



Conditions (i) and (ii) from Sec. [IJ arc thus fulfilled since 
the Fermi surface p 2 = Imp is fully gaped. 

When (#3 = 0, the auxiliary single-particle Hamiltonian 



( |3.13 ) represents two independent copies of the original 
chiral p-wave BdG Hamiltonian (3.10). When </> 3 = 0, the 
spectrum of Hi ux is identical to the spectrum of 7/? d 5„ 
up to a two-fold degeneracy arising from the doubling. 
Furthermore, <p 3 = implies that, in addition to the 
antiunitary charge conjugation symmetry 



c P "h «r(p> <^ x )) c Ph = - n T x (p> ( 3 - i5a ) 



with 
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where K denotes as before the operation of complex con- 
jugation, that originates from Eq. (3.12), there exists an 
auxiliary unitary charge-conjugation symmetry 

C-jH^ip, C ch = -H; ux (p, 4>(x)) (3.16a) 

with the generator of the auxiliary chiral transformation 



Po 
Po 



C ch := I" ™). (3.16b) 



Although neither C ph nor C ch are symmetries as soon as 
4> 3 7^ 0, their product 7~ aux = C ph C ch is a symmetry for 
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any <j> 3 . The antiunitary operation 7^ ux can be thought 
of as implementing an auxiliary time-reversal symmetry. 
As a result all eigenstates of %^ ux (p, <f)(x)), including 
zero modes, have a two-fo ld Kr amers degeneracy. 

As we shall see in Sec. [II C , Dirac Hamiltonians can 
also share unitary and antiunitary charge-c onjug ation 



symmetries. There are differenc es wi th Sec. |III Q how- 
ever. A first difference with Sec. Ill C is that 



C*H™(p, </>(*)) C ch = -H°T(P, C ch cf>(x)) (3.17a) 



with 
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A second difference with Sec. Ill C is the absence of rcl 



ativistic invariance. A third difference is that Eq. (|2.23| ) 
is non- vanishing. 



We start from the expansion (|2.1S| ) of the quasiparticle 
charge density induced by the order parameter, the static 
triplet (p. Wc compute first the contribution from n = 
1 for the pair of auxiliary Hamiltonians. It vanishes. 
The adiabatic approximation to the quasiparticle charge 
density for any one o f the pair of auxiliary Hamiltonians 
is given by Eq. fl2.22p 
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When fi < 0, Eq. ( p.18; ) gives p^*(x) = 0. This is con- 
sistent with the absence of a normalizable zero enersj^soj 
lution for negative values of the chemical potential.E2rc3 

When /i > and 3 is independent of space (a;), the 
adiabatic approximation to the conserved quasiparticle 
charge of the auxiliary Hamiltonian is 
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where the parametrization 



limit in which (f> 3 tends to zero, we get for the induced 
charge of the auxiliary Hamiltonian 



dia = ± ^ / 



(3.20) 



± winding number in (t^, </> 2 ). 



To compute the number A au x of u noccupied zero modes 
with the counting formula (2.12) induced by a vortex 
with vorticity unity, we choose the charge-conjugation- 
symmctry-brcaking perturbation such that it shifts the 
zero mode in energy above the chemical potential. In the 
limit 3 — y 0, we find that 



(3.21) 



Equation (3.21) implies that the adiabatic approximation 
for the number N of un occup ied zero modes of the orig- 
inal BdG Hamiltonian ( |3.10| ) induced by a vortex with 
unit vorticity is 



N= 1. 



Chern number 



(3.22) 



Next, we com pact ify space ( |3.1c ) and the order- 
parameter space ( |3.1a ), 
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where we have introduced the spherical coordinates = 
(6i,9 ? ) o f the two-sphere S 2 C R 3 . Motivated by 
Eq. ( 2.22 ), we separate the adiabatic approximation to 
the conserved quasiparticle charge of the auxiliary Hamil- 
tonian into two contributions, 



Qadfa = Q 
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(3.19b) 



is assumed for the superconducting order parameter with 
r and 6 denoting the polar coordinates of x G M. 2 . In the 



and g := (icj — %^ ux ) . Th e sub script refers to the 



semi-classical Green function ( 2.20 ). The first contribu- 
tion comes about when the order parameter decouples 
from the momentum. The second contribution arises 
when the order parameter and the momentum couple. 
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By explicit computation of the trace in Eq. (3.24b), 
one verifies that this trace is fully antisymmetric in all 
the indices i 2 , a 2 , «i, a x where the family i and the family 
a of indices are distinct with 



The integer 1? is the smallest dimensionality compatible 
with an irreducible representation of the Clifford algebra 
generated by the unit R x R matrix 1 and the traceless 
and Hermitian matrices 



1,2, 



1,2,3. 



(3.25) 



Furthermore, the trace in Eq. ( 3.24b| ) gives an integrand 
with even contributions in uj so that it does not vanish 
upon integrati on ov er ui. By explicit computation of the 
trace in Eq. ( [3.24c ), one verifies that this trace yields 
an integrand that is odd in u> and thus vanishes upon 
integration over oj, 



rt 



28,. „1, (3.28b) 
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where fj,,v = i,..-,2d + lE We choose the single- 
particle Dirac Hamiltonian such that a i?-dimcnsional 
spinor is coupled to a (d + l)-tuplet of real-valued static 
Higgs fields 
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S = o. 



(3.26) Accordingly, 



Wi th th e hel p of the manipulations made between 
Eqs. (C34) and ( C45 ), it is possible to write 
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where we have introduced the family of indices labeled 
by v 



1, 



,5, 



the five-momentum 
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and the domain of integration 
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(3.27d) 



The "surface" element of the sphere S 2 is here denoted by 
df2 2 - The subscript refers to the semi-classical Green 



function— (pT20[) . Equation (3.27) is the second Chern 
!r£3lt 



number.ni It takes integer values only. It should be com 
pared with Eq. ( p28Q 
andlcT 



We defer to Sees. QIC, [II D 



for a more detailed discussion of the connection 
between this second Chern number and the number of 
unoccupied zero modes. 



We have verified that the second Chern number (3.27) 
vanishes for a s-wave superconducting order parameter 
that supports an isolated vortex with unit vorticity. 



d 

H° irac (p,0(x)) :=^r^ 
1=1 



d+1 
a = l 



The decomposition (2.5) is generalized to 
<(>(x) = 4> + Sd>(x) 



[x). (3.28d) 



(3.29) 



and is assumed to hold together with conditions (iii) from 
Sec.0 once we have established the existence of charge- 
conjugation symmetry. Conditions (i) and (ii) follow 
from the spectrum 



4(p) 



(3.30) 



ofW° irac (p,0 o ). 

The Dirac Hamiltonian (3.28d) has the charge- 
conjugation symmetry 



-2,/Dirac 
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d+a 



'1/DiraCT^ 
n d 1 d+a 



(3.31a) 



for any a = l,---,d+las soon as (f> a = everywhere in 
space. Without loss of generality, we choos e the genera- 
tor of the so-called chiral symmetry ( 3.31a ) to be r 9f , +1 , 
i.e., we identify the charge-conjugation symmetry (2.4c) 
with the unitary operation represented by 



C, 



ch 



2d+ 



1 = (-i) d T 1 ---T 
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Equation ( |3.3l| ) follows from 2d + 1 being an odd inte- 
ger E3_With this choice, the vector-valued order parame- 
ter (U) is 



(3.32) 



i.e., 



Dirac single-particle Wj lmc for any d 



Our pur pose here is to apply the counting for- 
mula ( |2.12| ) to a Dirac single-particle Hamiltonian de- 
fined in the rf -dimensional Euclidean space with the coor- 
dinates (|2.l| ). We choose the dimensionality of the Dirac 
matrices to be 



D 



(3.33) 



It is shown in Appendix |^ that ther e exist s a symmet- 
ric matrix T d in the Clifford algebra ( [3. 28b ) that imple- 
ments the time-reversal symmetry 



1/Dirac 

H-d 



= r d - 1 nf 
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(3.34) 



When 



R 



(3.28a) 
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10 



the s ymmet ries (3.31a) and (3.34) of the Dirac Hamilto- 
nian (3.28d) can be combined into the antiunitary oper- 
ation 



C 



P h 



T d KC, 



cli 



(3.36) 



that implements a particle-hole transformation and un- 
der which the Dirac Hamiltonian ( |3.28d| ) is invariant, 



{T d T 2d+1 ) 'nf^* {T d T 2d+1 ). (3.37) 



We present our main result for Dirac fermions in three 
steps. The reader is referred to Appendix |c] for their 
detailed derivations when d = 1 and d = 2, and a sketch 
of the derivation when d > 2. 

First, the adiabatic approximation to the counting for- 
mula (2.12) is determined by the conserved quasiparticle 
adiabatic charge Q adia - However, <2 adia can be computed 
without assuming charge-conjugation s ymm etry. It is 
only when using the counting formula (2.12) that <3 adia 
must be restricted to a charge-conjugatio n-sym mctric 
configuration of the Higgs order parameter ( |3. la ). After 
relaxing the condition of charge-conjugation symmetry, 
we find 

Qadia = Chern d 

(-l)d(d-P/2(-i)*+i dl /o r 

(2d+l)\ ( ' J ^-^+1 

K 
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(3.38a) 



if we compactify momentum space and the target mani- 
fold, i.e., we have introduced the Euclidean momentum 



(3.38b) 



with (p 1: - ■ ■ ,p d ) £ S d and the spherical coordinates 6 
(<?!,••• ,9 d ) with 0(0) G S d C 



The last inclusion 
serves to emphasize that the coordinates on the d-sphere 
in order-parameter space can involve a sizable breaking 
of the conjugation symmetry through 4> d+1 (0) = (f)(0). 
The short-hand notation J stands for the integral over 

K 

R x S d x S d defined in Eq. ( |C45a| ). The Euclidean single- 
particle Green function g = (iui — H^ lrac ) 1 is the one 
for the Dirac Hamiltonian. T he su bscript refers to the 
semi-cla ssical Green function ( 2.20 ). The second equality 



in Eq. (3.38a) defines the d-th Chern number in term 
of these Dirac Euclidean single-particle Green functions. 
It takes integer values and, conversely, to each integer 
-^hedgehog there corresponds a static <p N (x) £ S d 

for which Q adia = 7V hcdgohog . We call 4>^2.(x) £ S d 
a hedgehog with the d-th Chern number N hedEehog . 

Second, if we rela x the condition that 4>{9) £ S d and 
replace it with Eq. (3.1b) instead, then the local quasi- 
particle current induced by the adiabatic variations of 



the Higgs fields in space (a;) and time (x ) is 



Q-l e U Vl -U d 



e^-a^A^ ■ • • d » d <t>* d - (3.39a) 



J ' \(p\ d ^ d\ 
It obeys the continuity equation 

d v f = 0. (3.39b) 

Here, 



101 := M + ~- + 



(3.39c) 



the (d + 1) indices v, v 11 ■ ■ ■ ,v d run over space (a;) and 
time (x Q ), and the (d + 1) indices a l7 --. , a d+1 run 
over the (d + 1) components of the vector-valued Higgs 
field ( |3.28c ). The Minkowski metric is used in space (x) 
and time (x ). Final ly, Q.^ denotes the area of the en- 
sphere S d . Equation ( 3.39a ) was obtained by Goldstone 
an Wilczck in Ref. |31| when d = 1 and d = 3 and by 
Jaroszewicz in Ref. |32| when d = 2. The conserved quasi- 
particle charge induced by a hedgehog configuration is 
the integer 



Qadia — J d d x j adia (a;) 



dO„ 



(3.40) 



s d 



■ "^hedgehog £ '^ J - 



The "surface" element ( "area" ) of the sphere S d is here 
denoted by 6D. d (£l d ). 

Third, we need to define a point defect that is compat- 
ible with the charge-conjugation symmetry and consis- 
tent with the adiabatic approximation (see Sec. Ill D for 
a discussion of the latter caveat). A charge-conjugation 
symmetric point defe ct is a half hedgehog, i.e ., a static 
Higgs (d + l)-tuplet (3.1a) that satisfies Eq. (3.1b) and 
wraps the d-sphere S d the half integer ±1/2 number of 
times (a domain wall in d = 1, a vortex wit h unit vor- 
ticity in d = 2, etc). The counting formula ( 2.1 2|) pre- 
dicts that there is one unoccupied zero modes for a half- 
hedgehog provided we choose the sign of the infinitesimal 
conjugation-symmetry-breaking <fr in Eq. (3.1a) to be op- 
posite to the sign of ±1/2. 

The d-th Chern number (3.38) is non- vanishing for 



any single-particle Hamiltonian topologically equivalent 
to the Dirac Hamiltonian ( [3. 28] ). 



D. Physical interpretation of the adiabatic 
approximation 

The relevant length scales in <i-dimensional space are: 
(1) the characteristic linear size a of a point defect, (2) 
the characteristic linear size £ = 1/A of the support in 
space of the zero modes bound to the point defect, (3) 
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FIG. 2: (Color online ) T he dependence on space (a;) of the 
Higgs doublet in Eq. (CI) is depicted in panels (a), (b), and 
(c) together with the number and position of the zero modes 
along the single-particle energy eigenvalue axis e. The mag- 
nitude of the asymptotic value of the order parameter cj) 1 
is ip > 0. The magnitude of the asymptotic value of the 
charge-conjugation symmetry (f> 2 is <f> > 0. The size of the 
point defect, a domain wall, is a. The exponential decay of 
the envelope of the bound state is controlled by the length 
scale £ n — 1/|</> |. The characteristic length £ is defined in 
Eq. full). 



the characteristic length £ over which the variation of the 
order parameter <fi is of order of the gap 2A , 



rw=A , 



(3.41) 



(4) and the linear size R of the region of space in which 
the conserved quasiparticle charge Q and the total num- 



ber TV of unoccupied zero modes is to be computed (mea- 
sured). The gradient approximation relics on the hierar- 
chy 



a < £ <C £ <C R. 



(3.42) 



In the gradient approximation, all microscopic data for 
length scales smaller than £ are dispensed with. This 
fact dictates how to properly interpret the results from 
the gradient approximation when point defects can be 
assigned an additive label q def . Here, the subscript stands 
for point defect. 

When d = 1, q dcf = q dw , with q dw = +1 representing 
a domain wall and q dc{ = — 1 representing an antidomain 
wall. Observe that we can exchange the terminology do- 
main and antidomain wall, for there is no absolute notion 
of a positive or negative q dw . When d = 2, q dci = q vm . 6 Z 
encodes the vorticity of the vortex. In either case, these 
charges can be thought of as "classical Coulomb charges" . 

Any two point defects with the labels q d J f and q d 2 J { 
within a distance of order a of each other fuse into a 
point defect with the label 



(1+2) _ (1) (2) 
?def ~ <?dcf ?def 



(3.43) 



When d = 1, domain walls are necessarily ordered. Con- 
secutive domain walls necessarily carry the opposite la- 
bel. Fusion of two domain walls in d — 1 necessarily 
results in their annihilation. When d = 2, vortices can 
fuse to yield larger vorticitics. 

When d = 1, we have shown in Sec. CI that the num- 
ber of unoccupied zero modes TV bound to a single do- 
main wall in a large region of size R is TV = 1 in the 
adiabatic approximation. This result agrees with solv- 
ing the differential equation for the gero mode once all 
microscopic data h ave b een supplicd.El However, we have 
also shown in Sec. [II A and Sec. C 1 that the number of 



unoccupied zero modes TV in a large region of size R can 
be made to be an arbitrary integer number in the adia- 
batic approximation, in apparent contradiction with the 
fact that the net number of domain walls is either or 
1 along any finite segment of the line. This is explained 
with the help of Fig. |^ as follows. 

Occupying or leaving empty any single-particle level 
is a physical operation forcing a fcrmionic quasiparticle 
in or out of this single-particle level. This physical pro- 
cess is, for zero modes, achieved by the local sign of the 
charge-conjugation-symmetry-breaking <f> 2 . In Fig. || (a), 
we plot as a function of position x along the line the 
dependence of the order parameter <p 1 and of the charge- 
conjugation-symmctry-brcaking <f> 2 - Given the shape of 
the domain wall and its asymptotic value ip , we choose 
the sign (f> 2 (x) = +4> so that the zero mode bound to the 
domain wall is shifted to a positive energy. If the chemi- 
cal potential is chosen to be at 0, then the zero mode is 
unoccupied. The zero mode bound to an antidomain wall 
remains unoccupied if we choose the sign 4> 2 ( x ) = ~ as 
depicted in Fig. |^(b). In Fig. ^(c) a domain wall is fol- 
lowed by an antidomain wall a distance £ apart. The clos- 
ing of the gap |0| at the center of each domain walls when 
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4> 1 = is avoided by the charge-conjugation-symmetry- 
breaking field saturating to its asymptotic magnitude 4>. 
A phase twist in cj> by w half way between the two con- 
secutive domain walls insures that the zero-modes are 
shifted to positive energies and thus remain unoccupied. 
Gap closing at this phase twist is again avoided because 
4> 1 has healed to its asymptotic value. Everywhere along 
the line, |V0| can thus be made as small as needed. The 
energy splitting 5e, which is of order A exp (— 2g/£ Q ), is 
exponentially small in £/£ . In the gradient approxima- 
tion, the order of limit is i/£ — > oo first followed by 
<f> — > 0. Evidently, this order of limit does not commute 
with 4> — > first followed by £/£ Q — > oo. In the latter 
order of limits, the unoccupied number of zero modes 
always vanishes. 

We close the discussion of the results 



m one- 



dimensional space from Sec. [II A by observing that since 



the number of unoccupied zero mode is an integer, all 
higher contributions to the gradient expansion beyond 
the adiabatic (leading) order must vanish identically. We 
have verified this expectation explicitly for the first sub- 
leading order. 

When d = 2, the adiabatic approximation of Sec. [II B 
applied to a chiral p-superconductor predicts that the 
total number of unoccupied zero modes in a large region 
of size R equals in magnitude the net vorticity in this 
region. On the other hand, it is known from solving the 
differential equation for a single vortex of vorticity <7 vort in 
a p-wave superconductor that the number xif zero modes 
is one if q VOTt is odd or zero otherwise. EjE3 This is not 
a paradox if the adiabatic approximation in this paper 
is limited to point defects each of which bind at most a 
single zero mode. We now argue that this interpretation 
of the adiabatic approximation is a necessary one. 

The adiabatic approximation is not sensitive to the 
microscopic data on length scales smaller than I. How- 
ever, it is precisely those microscopic data that deter- 
mines if more than one zero modes can be bound to a 
point defect. Consider the case of the relativisti c Di rac 
Hamiltonian in two-dimensional space from Sec. C2. It 
respects two charge- conjugation symmetries, the chiral 
symmetry and the particle-hole symmetry. Correspond- 
ingly, the order parameter ip .can cither be interpreted as 
a bond-order (Kekule order) Jj in which case the electron 
charge is a good quantum number, or as an s-wave su- 
perconductor^ in which case the electron charge is not 
anymore a good quantum number (the thermal quasipar- 
ticle charge is). In either interpretations, an index the- 
orem guarantees that the number of zero modes e quals , 
in magnitude the vorticity of the order parameter .P 4 lr 5 I 
For example, a vorticity of two implies that there are 
two zero modes. Of course, they must be orthogonal. 
Perturb now the Dirac Hamiltonian with a perturbation, 
whose characteristic energy scale is rj, that breaks the 
chical symmetry but preserves the particle-hole symme- 
try.E3 We assume that r\ is not sufficiently large to close 
the gap. The delicate balance that allowed the two zero 
modes to be orthogonal is destroyed by the perturbation 



rj. The two zero modes split pairwise, one migrating to 
positive energy, the other migrating to negative energy. 
This level repulsion is encoded by the microscopic data 
in a region centered about the point defect of linear size 
£ , a window of length scales inaccessible to the adiabatic 
approximation. If we split the single vortex with vorticity 
two into two vortices with vorticity one separated by a 
distance £, we can use the adiabatic approximation. The 
energy shift induced by the chiral-symmetry-brcaking <p 
is of order </>exp(— 2£/£ ). The level splitting induced by 
the particle-hole-symmetric rj is of order ?yexp(— 2£/£ ). 
They vanish in the adiabatic limit £/£q — > oo first, — > 
second. 

The adiabatic approximation saturates the number of 
unoccupied zero modes N. Consistency demands that 
higher-order corrections vanish identically. This suggests 
that the adiabatic approximation cannot capture the par- 
ity effect by which a charge-conjugation-symmetric per- 
turbation splits pairwise the degeneracy of zero modes. 
Such a parity effect is an essential singularity for the adi- 
abatic expansion presented in this paper. 



IV. 



CONCLUSION 



In this paper, we have derived a procedure to count the 
zero-energy eigenvalues of a single-particle Hamiltonian 
H(p,(p(x)) that possesses a charge-conjugation symme- 
try, when the order parameter tp(x) supports point de- 
fects. This procedure is based on counting the conserved 
quasiparticle charge Q induced by a point defect of <p(x). 

The conserved quasiparticle charge Q exists because of 
the Hcrmiticity o{H(p,ip(x)y For instance, in problems 
for which the electric charge is a good quantum number, 
the number of zero modes can be determined from the 
induced electric charge near point defects. In supercon- 
ductors for which the mean-field approximation is accu- 
rate, the conserved charge is associated to the thermal 
currents carried by the Bogoliubov quasiparticles. 

Irrespective of the origin of the conserved quasiparticle 
charge (electric or thermal) , we showed that one can use 
counting arguments similar to those appearing implicitly 
in Ref. ^|and explicitly in Ref. ^to compute the conserved 
fractional quasiparticle charge Q induced by a point de- 
fect. Once one observes that the number of unoccupied 
zero modes N can be related to Q through N = —2Q, 
one can concentrate the efforts into computing Q near a 
point defect of the position-dependent order parameter 
<p(x). We carry out this procedure within a gradient ex- 
pansion for smoothly spatially varying fluctuations of the 
chargc-conjugation-symmctric order parameter (p(x)E3 

The resulting expression is then applied to generic 
systems in one-dimensional space, the chiral p-wave su- 
perconductor in two-dimensional space, and to Dirac 
fermions in d-dimcnsional space. In one-dimensional 
space, we find that the number of unoccupied zero modes 
is related to the first Chern number. For the p-wave su- 
perconductor in two-dimensional space, it is related to 
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the second Chern number. For the Dirac Hamiltonians 
in cf-dimensional space, the number of zero modes is de- 
termined by the d-th Chern number. Therefore, we can 
establish in a logical and constructive way the relation 
between the number of zero modes induced by a point 
defect and topological invariants (the Chern numbers) in 
a number of cases. 



There has been a resurgence of efforts dedicated to 
counting zero modes induced by defectiv e ord er parame- 
ters in the condensed matter community.! 36 !! 37 ! Fukui and 
Fujiwara in Ref. |3^ have revisited Dirac fermions in up 
to three spatial dimensions. They use the chiral anomaly 
to carry the counting. Our general counting formula re- 
produces their and previous results, and extends them to 
arbitrary d dimensions of space. An elegant formulation 
of the Dirac problem in d dimensions in the presence nf. 
an isotropic hedgehog has been carried out by Hcrbut.EHl 
(See also Freedman et al. in Ref. ^0[) The number of zero 
modes should not depend on deformations away from the 
isotropic case, and this result is captured by our counting 
formula, which we can express as the d-th Chern num- 
ber. Finally, Teo and Kane in Ref. |l] have studied de- 
fects of arbitrary dimensions r coupled to nonintcracting 
fermionic quasiparticles with the help of the classifica- 
tion scheme of Schnydcr ei_qLfbr topological band insula- 
tors and superconductors .CjCj They conjecture that the 
topological invariant for a given symmetry class is related 
to the number of zero modes attached to r-dimcnsional 
defects, but cannot make any direct connection between 
these two integer numbers. Our explicit construction pro- 
vides this relation, although, only for specific examples 
and for point defects. We do not have yet a proof for 
generic Hamiltonians that have a momentum-dependent 
coupling to the order parameter. 



In summary, we showed that one can count zero modes 
induced by point defects using a conserved quasi-particle 
charge present for any Hermitian mean-field Hamiltonian 
with charge- conjugation symmetry. This counting sup- 
ports a direct relation between topological invariants and 
the number of zero modes bound to defective order pa- 
rameters. 
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Appendix A: Counting the zero modes with the 
conserved quasiparticle charge 

In this Appendix, wc review an identity that relates 
the total number of zero modes of any single-particle 
Hamiltonian % with charge-conjugation symmetry to its 
Euclidean Green function Q (u>) := (iuj — H)^ 1 . This ider 
tity appears implicitly in Ref. ^ and explicitly in Ref. 
We will only assume that the single-particle Hamiltonian 
T-L has the property that there exists a transformation C 
such that 



c^nc = -U, 



(Al) 



whereby the transformation C is norm preserving, i.e., it 
can cither be an unitary or an antiunitary transforma- 
tion. We call C the operation of charge conjugation. 
To simplify notation, we take the spectrum 



{0, sgn(n)e H | n = ±1, ±2, • • • } 



(A2) 

of H to be discrete up to finite degeneracies of the eigen- 
values. This is the case if H describes the single-particle 
physics of a lattice Hamiltonian. The spectral decompo- 
sition of TL thus reads 



«= E 



> o, 



S S n ("-) £ |n|l ? /'n}( ? /'n|, 
n=±l,±2,--- 

with the single-particle orthonormal basis obeying 



(A3) 



(4>n\4>r. 



N 



Ei«)( c 



(lpn\oc) = 0, 
E \^n}{lpr 



(a\a') 



(A4) 



n=±l,±2, 



We have assumed the existence of N zero modes labeled 
by the index a. The relation 



\ip n ) = C\ip- n ) 



(A5) 



holds for any finite energy eigenstate labeled by the in- 
dex n = ±1, ±2, ■ • • as a result of the charge-conjugation 



symmetry (Al). 

On the lattice, wc denote the value of the energy eigen- 
functions at site i by 

V' sg n(n) E|n|ii : = (#n), n = ±1, ±2, • • • , (A6) 

for the finite-energy eigenvalues and 

1> ati :=(i\a), a =1, (A7) 
for the N zero modes. For any two lattice sites i and j, 

= E ^-i„i,^_i„i,i 

n=-l,-2,— 

+ E ( C ^-\n\,j (WUn\,i 
n=-l,-2,--' 

N 



(A8) 



E^«.j^W' 
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wher e wc have used the completeness relation defined in 
Eq. (A4). When i = j, we find the local sum rule 



N 



1 = 2 J2 r-\ n \,^-\n\,i + J2<^i (A9) 



owing to the fact that C is norm preserving. 

We now assume that there exists a single-particle 
Hamiltonian 



n=±l,±2,-- 

= — C ^ 1~Lq c 



(A10) 



such that H does not support zero modes. Evidently, the 
local sum rule 

1 = 2 ]T 4>-\n\,^-\n\,i (AH) 
n=-l,-2,--. 

also applies. 

After subtracting Eq. ( All ) from Eq. (A9), the local 
local sum rule 

-2 (V'-|„|,iV'-|n| 1 » -^-|n|,t^-|n|,i) = 

n=— 1,— 2,— 

n (A12) 

a = l 

follows. In turn, after summing over all lattice sites and 
making use of the fact that the zero modes are normalized 
to one, the global sum rule 



iV 

T 



-|n|,t^-|n|,i ^-H.^-Mi 

n=-l,-2,— i 

(A13) 

follows. The global sum rule is only meaningful in the 
thermodynamic limit after this subtraction procedure has 
been taken. 



We now present the global sum rule ( A13 ) with the 
help of Euclidean single-particle Green functions. Wc 
thus define the Euclidean single-particle Green functions 



and 



1 



iu-H 



(Al4a) 



(A14b) 



for any real- valued and non- vanishing uj. Next, we choose 
7 to be the contour in the complex uj plane that runs 
counterclockwise along the real axis ui G M, avoids the 
origin uj = by an infinitesimal deformation into the 
upper complex plane Re a; > 0, and closes through a 
semi-circle in the very same upper complex plane (see 



CO complex plane 




FIG. 3: (Color online) Definition of the integration contour 
7 that picks up the discrete negative energy eigenvalues of 
the single-particle Hamiltonian. The shaded box represents 
the Fermi sea. The filled circles represent the discrete energy 
eigenvalues. The contour 7 avoids the zero mode. 



Fig. g). With the help of the residue theorem, it then 
follows that 



n=— 1,-2,— 



dui 
2^ 



[i\g(u)\i) (Al5a) 



and 



n=— 1,— 2,— 



did 

2^ 



■i\G (u)\i). (A15b) 



Equation ( A13) can now be rewritten in the desired form 



Y = ~ Tr / tjrIW-SoH]- (Al6a) 



Equation (Al6a) can be expressed in terms of the local 
"lattice charge" p i through 



N 



2^ 



(i\[g(u)-Q (u)]\i), 



(A16b) 



or the local "continuum charge" through "charge" p(r) 

f = -/ d d rp(r)^-Q, 

f dui (Al6c) 
Pir) := / 2^(r|[0M-S o ( W )]|r>. 
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Hcrmiticity of T-L and H guarantees that these "quasi- 
particle charge densities" obey a continuity equation, i.e., 
that their quasiparticle charge Q is conserved. 



Appendix B: Time-reversal symmetry of 'H^ >lrac 



We are going to prove the existence of a time- 
rever sal s ymmetry for the Dirac Hamiltonian defined by 
Eq. (3.28). To this end, we shall introduce an auxiliary 
Lagrangian £ Dirac ift,(2d+l)-dimensional space and time. 

By construction,!^ 



£-2d+i( X l' " ' > x 2d' X 2d+l) 



2d+l „ 

£ r «£ (B1) 

u— 1 ,y 

can be interpreted as the free Euclidean Dirac Lagrangian 
in (2d + l)-dimensional Euclidean space (x) and time 
( x 2d+i)- Observe here that 2 d , where d = [(2c? + l)/2] 
is the smallest integer larger or equal to the real number 
(2d+l)/2, is the smallest dimensionality of an irreducible 
representation of a Dirac Hamiltonian in 2d-dimcnsional 
space.123 In the reminder of the argument, we adopt the 
summation convention over repeated indices. 
The classical action 



.5' 



xip a (x) r 



(B2) 



v / a/3 



is invariant under the group 0(2d + 1) generated by 



1 < v < y' < 2d+ 1. 



(B3) 



Consequently, there must exist (2d 4- 1) unitary matrices 
V v from the (complex) Clifford algebra such that 

(B4) 

Symmetry under reflection about any direction v = 
1, ■ • • , 2d+ 1 in (2<i+ l)-dimcnsional space and time thus 
becomes 



/-Dirac/ 

-1 /-Dirac/ i 



2d+l 



(B5) 



There must also exist an auxiliary unitary matrix T av 
from the (complex) Clifford algebra under which 



£-2d+l( X li ' ' ' i X 2d' X 2d+l) 
Taux -^Dirac (~^~ x l > ' ' ' 



b 2di 



L 2d+l) Taux 



We are now ready to define the unitary matrix 

d 

7~d '■= T aw J > 2d+l Y\. ^ a 



(B6) 



(B7) 



from the (complex) Clifford algebra. If K denotes com- 
plex conjugation, we interpret the operation T d K as 
imple mentin g reversal of time on the Dirac Hamilto- 
nian (3.28d). It is a symmetry, for 

H»™ = T d - 1 H»™*T d . (B8) 



Appendix C: Zero modes of irac 



T his A ppendix provides intermediary steps for 



Sec. IIIC 



Dirac fermions in one-dimensional space 



a. Definition 



When d = 1, Eq. (3.28) simplifies to 



H^{p,<j>(x)) -a^p + a^^+a^x). (CI) 

The (R = 2)-dimensional representation of the Clifford 
algebra is here generated from the Pauli matrices a 1 , a 2l 
and a 3 . If the doublet of Higgs fields <p Q is constant 
through one-dimensional s pac e (x), the single-particle 
spectrum of Hamiltonian (CI) has a gap controlled by 
the Higgs components adding in quadrature, 



4(P) = P 2 + 00 = P 2 + 4% 1 + 00 2- 



(C2) 



The generator of the chiral symmetry of the Dirac 
Hamiltonian (CI) is 



if 



C ch = °2 



02 =0 



(C3a) 



(C3b) 



everywhere in Euclidean space i£l 

The operation of time-reversal is implemented by 

T x K := <Tj K (C4) 

where K denotes complex conjugation. It is a symmetry 
of the Dirac Hamiltonian (CI) for any Higgs configura- 
tion and <p 2 ■ 

The operation of particle- hole exchange is implemented 



by 



Cph^TiKC, 



ch 



-icr 3 K 



(C5) 



and it is only a symmetry of the Dirac Hamiltonian ( |Cl| ) 
provided (j) 2 = 0. 

The discovery that this model supports zero modes 
was made by Jackiw and Rebbi in Ref. [&]. Its relevance 
to the physics of polyacetylenc was made by Su, Schri- 
effer, and Heeger in Refs. || and [?]. In polyacetylene, 
the kinetic energy results from linearizing the dispersion 
around the two Fermi points of a single-band nearest- 
neighbor tight-binding model at half-filling in the left- 
and right-mover basis. The Higgs field (f> 1 = (p 1 real- 
izes a modulation of the nearest-neighbor hopping am- 
plitude that is mediated by phonons. A Peierls tran- 
sition opens up a single-particle electronic gap through 
the breaking of the translation symmetry by one lattice 
spacing down to a residual translation symmetry by two 
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lattice spacings. A domain wall in tp 1 that interpolates 
between the two possible dimer ground states binds one 
zero mode per electronic spin (which we have ignored so 
far). The Higgs field <p 2 breaks the sublattice symme- 
try of the tight-binding model. The quasiparticle charge 
density induced by a domain wall with a non-vanishing 
charge-conjugation symmetry breaking <j) 2 was computed 
by Goldstone and Wilczek in Ref. ||l] and shown to vary 
continuously with <f> 2 . We are going to repro duce all these 
results using the adiabatic approximation (2.21). 



Counting zero modes 



We start from the expansion ( 2.18a ) of the quasipar- 
ticle charge density which, for a Dirac Hamiltonian, is 
exact and consider the contribution from n = 1. The 
adiabatic approximation to the quasiparticle charge den- 
sity is 



Padia( x ) 



dudp 2e ah {d x 4> a ) {x) (j) h {x) 
(d x </> a )(x)cf> b (x) 



2tt [#(z) + 



For a constant <fr 2 > that breaks the conjugation sym- 
metry, integration over one-dimensional Euclidean space 
(x) gives the adiabatic approximation to the conserved 
quasiparticle charge 



b 1 (x— +oo) 



fadiaVW 



2tt {eft + <ft) 



X — — OQ ) 



(C7) 



(27r) 1 ( arctan 



x—-\-oo 



The single domain wall with the asymptotic values 

<M±°°) = = FPo> ¥>0> : ( C8 ) 

here chosen in such a way that any zero mode is shifted 
in energy by <f> 2 > above the chemical potential, induces 
the negative conserved quasiparticle charge 



l im + <3adia<>2) = 

o ^o+ 2 



(C9) 



Having restored th e ch arge-conjugation symm etry b y re- 
moving cj) 2 in Eq. (|C9|), the counting formula (2.12) can, 
in turn, be used to deliver the positive number of unoc- 
cupied zero modes 



N = 1 

bound to this single domain wall. 



(CIO) 



Chern number 



Whereas the counting formula ( 2.12| ) relies on the 
charge-conjugation symmetry, the adiabatic approxima- 
tion to the quasiparticle charge density does not. We 
are going to take advantage of this freedom to relate the 
conserved quasiparticle charge to the first Chern number. 

To this end, we compactify the base space x G R to the 
circle x € S 1 by imposing periodic boundary conditions. 
We then parametrize the doublet of Higgs fields according 
to the polar decomposition 



<P(x) 



eb 2 {x) 



sin 0{x) 
cosO{x) 



(Cll) 



We have thus compactified the target space <p(x) 6 M 2 
to the unit circle 9(x) S S 1 . 



In the adiabatic approximation (2.21), the conserved 
quasiparticle charge is 



1 



* adia 



2(2tt) 



da; 



2T! 27T 



x tr g^—g^L—g - gJL—Q-t—g 



dp 86 



86 dp 



(C12) 



The integrand can be written in a more compact and 
symmetric form by observing that 



dg- 1 

du> 



o 



. dg- 1 

duj 



o \ / o 

and by introducing the three-momentum 

K v = (u,p,6) e R x S 1 x S 1 . 



(C13) 



(C14) 



Equation ( C12 ) is then nothing but the first Chern num- 
&er,E3 for 



-1 



^adia " 



d 3 Ke M „ A tr 



dK„ dK„ dK\ 



24tt 2 , 

KxS 1 xS 1 

_ (C15) 

We infer that the conserved quasiparticle charge induced 
by the adiabatic winding of the Higgs doublet around the 
circle takes integer values. 



Moreover, the domain wall from Sec. |C l| is a half- 
winding of the unit circle S . More precisely, evalua- 
tion of the trace in the integrand of the first Chern num- 
ber (CIS) gives the adiabatic approximation to the con- 
served quasiparticle charge 

2tt 



i adia 



= winding number in <$. 



(C16) 



The integrand is nothing but the space (a;) and time (t) 
dependent quasiparticle charge density 

1 



J 



2tt|0| 



2 abra 



d,. 



(C17a) 
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that obeys the continuity cquationEl 

dS = o 

with 

1 .. 



2vr0 : 



,6,d,,,(h,,. 



2 aa Va^i/' Va 



(C17b) 



(C17c) 



Dirac fermions in two-dimensional space 

a. Definition 



When d = 2, Eq. Q3.28| ) simplifies to 

HY" ac {p, <t>(x)) := a^j + a 2 p 2 + + /3 2 <£ 2 + /3 3 3 . 

(Cl8a) 

The (i? = 4)-dimensional representation of the Clifford 
algebra can be chosen to be generated from the five trace- 
less and Hcrmitian matrices 



r, 



(01)02)^1)^2)^3) 



with 

and 

Pi ■■= *1 



: = a 2 



'2) 



(C18b) 
(C18c) 



(C18d) 

A second set of Pauli matrices 7\, r 2 , and t y has been 
introduced, together with the unit 2x2 matrices a and 
t . If the triplet of Higgs field <f> is constant throughout 
two-dimensional spac e (a;), the single-particle spectrum 
of Hamiltonian (C18) has a gap controlled by the Higgs 
components adding in quadrature, 



4(p) =p 2 + 4>1 = p 2 



*§1 



*0 2 



Aga- 



(C19) 



The gener ator of the chiral symmetry of the Dirac 
Hamiltonian (CIS) is 



if 



C ch = & = °3 







(C20a) 



(C20b) 



everywhere in Euclidean space x G M 2 . 

The operation of time-reversal is implemented by 



T-2 K := o~ x ®t x K 



(C21) 



where K denotes complex conju gation. It is a symmetry 
of the Dirac Hamiltonian (C18) for any Higgs configura- 
tion <j> 1 , (f>2, and <f> 3 . 

The operation of particle-hole exchange is implemented 

by 



C 



P h 



r 2 Kc ch 



7 2 (g) T 2 K 



(C22) 



and it is only a symmetry of the Dirac Hamiltonian (CIS ) 
provided 4> 3 = 0. 

The discovery that this model supports zero modes 
(Majorana fermions) was made by Jackiw and Ro ssi in 
Ref. |3^ within an interpretation of Hamiltonian (CIS) 
as a relativistic superconductor. Weinberg shortly there- 
after proved an index theorem in Ref. p3 for these zero 
modes. The effect on the induced conserved quasipar- 
ticle charge by a triplet of Higgs fields was investigated 
by Jaroszewicz in Ref. [32] (see also Refs. [45|-fl7|). It was 
proposed by Hou, Chamon, and Mudry in Ref. |^ that 
graphene could realize Hamiltonian (C18) with the Higgs 
doublet 4> 1 and <fi 2 responsible for a Kekule bond-density- 
wave instability and the charge-conjugation-symmetry- 
breaking (f> 3 responsible for a charge-density-wave insta- 
bility (see also Refs. 



Counting zero modes 



We start from the expansion (2.18a) of the quasipar- 
ticlc charge density induced by a static triplet of Higgs 
fields <p which, for a Dirac Hamiltonian, is exact. We 
compute first the contribution from n — 1. It vanishes. 
The adiabatic approximation t o the quasiparticle charge 
density is in fact given by Eq. ( 2.22| ) 



dujd 2 p 8e abc (d^J (x) (<9 2 b ) (x)(/) c (x) 



(2tt)= 



■ p 2 + (I> 2 (x)Y 



r abc 



(e^J (x) (d 2 4> h ) (x^x) 



4- 



\<P(x)\ 



(C23) 



For a constant <p 3 , the integration over x gives the 
adiabatic approximation to the conserved quasiparticle 
charge 



"0 

QadiaW> 3 )=^ / de J dpp 



(p 2 + tif 2 



where the parametrization 



p(r) cos6(6>)\ 
p(r) sin 9(0) J 



(C24a) 



(C24b) 



is assumed for the charge-conjugation-symmetric doublet 
of Higgs fields with r and denoting the polar coordi- 
nates of x € 
the induced charge 



In the limit in which (f> 3 tends to zero, 



' adi; 



») = 



sgn(0 3 ) 



x winding number in ((f> 1 , <p 2 ) 

(C25) 

follows. To compute the number of unoc cupied zero 
modes N with the counting formula ( [2.12 ), we choose 
the sign of the charge-conjugation-symmctry-brcaking <j) 3 
such that it shifts the zero mode in energy above the 
chemical potential, i.e., with the opposite sign to the 
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winding number of the Higgs doublet ((/) 1 ,^ 2 )- We then 
take the limit (f> 3 — > 0. If so, for a unit winding 



N = 1. 



(C26) 



Observe that the number of zero modcsJC24d) agrees 



sJC24a 
Ipppl 



with Weinberg's index theorem in Ref. |3J applied to a 
single vortex with unit vorticity. 



c. Chern number 

We use the notation Q := (iuj — H®"^)" 1 an d com- 
pactify both space and the order-parameter space, 



x £ S 



4>{0) e s 2 c 



(C27) 



where 6 = (0 1 ,9 2 ) are the spherical coordinates on the 
two-sphere S 2 cM 3 . 



With the general manipulations of Appendix G 2 , it is 
then possible to write 



Q 



adi, 



-K27T) 

60 



tr 4 (Qd, J g- x ---gd, J _g- 1 



K 



where we have introduced the family of indices v 



the momentum 

K v = (^Pi^AA) 
and the domain of integration 

J = J— J dn2{P) f dn2{6) 

K 



(C28a) 
(C28b) 
(C28c) 



(2tt)2 J (2tt)2 ■ (C2M) 



pes 2 



ees 2 



The "surface" element of the sphere S 2 is here denoted by 
df2 2 - The subs cript refers to th e semi-classical Green 
function^ ( |2-20| ) . Equation (C28) is the second Chern 
numberSB It takes integer values only. 



Chern number for Dirac fermions in 
d-dimensional space 



To prove Eq. (3.3S), imagine that we integrate out the 
Dirac fermions in the ba ckgrou nd, not necessarily static, 
of the Higgs fields in Eq. ( 3.28d ) subject to the constraint 



0(0) G S d C 



pd+1 



(C29a) 



where 9 are the polar coordinates of the d-sphere S d C 
The conserved current of the fermionic single- 
particle Dirac Hamiltonian must induce a conserved cur- 
rent Jadia w hh v = 0, 1, • ■ • ,d for the Higgs fields. Its 



time-li ke co mponent j^dia enters in the counting for- 
mula (2.12). Relativistic covariance, current conserva- 
tion, dimensional analysis, and the constraint (C29a) all 
conspire to bring this current to the form 

Here, summation convention over the repeated indices 



! "d 



0,1,--- ,d, 



• • • , a d, a d+i = d + i2d+ 1, 



(C29c) 



is implied. (Compared to our convention in the defini- 
tion (3.28d) of the Dirac Hamiltonian, we have shifted the 
values taken by the family of indices a = d+l,---,2d+l 
to stress that it differs from the family of indices i = 
,d.) 



If we co mpare Eq. (C29) with the gradient expan- 
sion ( 2.24 ), we deduce that the leadi ng no n-vanishing 
contribution to the gradient expansion (2.24) must be of 



order n = d and, for a static Higgs background, given by 

= H)%a,.. hai (\K) ■ ■ ■ (*<A) («)■ 

(C30a) 

The summation convention over the two distinct families 
of indices 



' ' ' > *d ~~ 1) ' ' ' l^l 

a l7 ■ • • , a rf = d+ 1, • • ■ , 2d + 1, 
is implied and the expansion coefficients are 



(C30b) 



^d a d"'H a l ' * 



2 71 



tr. 



dPi, dcj) 3j 



duj 



(w,p). 
5 (C30c) 



Any permutation of the indices on the left-hand side is 
defined by the same permutation of the differentials in 
the trace of the right-hand side. 
We are first going to prove that 



Qadia '■ ~ J ^ X jadia( X ) 



(_)d(d-l)/2^_^d 



d\ 



d"d °d"d "1"1 



(0) 



9<£S d 



(C31a) 



where the target space S d C M. d+1 of the order-parameter 
is parametrized by the (i-independent spherical coordi- 
nates 



8 



y d+l' 



(C31b) 
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the integral 

J =Jdn d (C31c) 
ees d s d 

with dSl d the "surface" element of the d-sphere, and 



Ji 



/ dg- 1 dg- 1 

"\ d Pu y 89. 




(w,p). 
(C31d) 



Again any permutation of the indices on the left-hand 
side is defined by the same permutation of the differen- 
tials in the trace of the right-hand side. 



Proof. To evaluate the trace in the integrand ( C30c ) , the 
semi-classical Green functions are re-massaged so as to 
bring all the Dirac matrices in the numerator, 



1 



Ld 2 + p 2 + <p 2 (x) 



(C32) 



Multiplying out al l Gree n functions in the trace from the 
integrand in Eq. (C30c) yields in the numerator terms 
made of the product from 2c? T-matrices, 2d + 1 T- 
matrices, 2d+2j T-matrices, 2d+2j + l T-matrices, 
2d+2d T-matrices, and 2d+2d+l T-matrices. Any trace 
over an even number 2d + 2j of T-matrices is odd under 
uj — > —us since it comes multiplied by the power uj 2d + l ~ 2 i 
in the numerator. Such a trace does not contribute to the 
uj integration since the denominator is an even function 
of u>. Any trace over an odd number 2c? + 2j + 1 of T- 
matrices is even under u> —lu since it comes multiplied 
by the power w +1 ~ 2 i~ l . Such a trace can only be non- 
vanishing if • • ,r 2d _|_ 1 _all appear in the trace and all 
an odd number of times.c3 For such traces, the Clifford 
algebra delivers another key identity in that 



tr, 



(C33) 



if u. = 1, • • • , 2c? for j = 1, • ■ • , 2c!. Here, R = 2 d and 



r 2d+1 was defined in Eq. ( 3.311: ). The c oefficients ( C30c ) 
inherit the antisymmetry of Eq. (C33) in that they are 
f ully a ntisymmetric under any exchange of the indices 
(|C30bD . 

We need to overcome the fact that the ranges of i and 
a are unequal in cardinality. To this end, we change 
variables on the target space and introduce the spheri- 
cal coordinates (C31b) of the targe t space S d C K d+1 . 
The adiabatic approximation ( p30a ) to the quasiparticle 
charge density becomes 



8; 



(*). 

(C34) 



(C35) 



Now, it is the summation convention over the two distinct 
families of indices 

*ii ' ' ' i id = 1j ' ' ' ' ' 
bi, • • • ,b d = c? + 1, - • • ,2c!, 

that rep laces ( C30b| ), whereby the expansion coef- 
ficients (C34) are related to the expansion coeffi- 
cients ( C30c| ) through the chain rule for differentiation, 
i.e., 



4 b 



d u d t l u l 



(C36) 



Here, any permutation of the indices of J on the left- 
hand side is defined by the same permutation on the in- 
dices of I on the right-hand side. 

By linearity, the antisymmetry ( C33| ) thus carries over 

to 



Ji 



'd"d '1"1 



(C37) 



MM 



(C38) 



where J\f d is a normalization constant. We will not need 
the explicit dependence of the normalization N d . We will 
only make use of the fact that it obeys the identity 

(dj)2 6 » d h d ■ ■ ■ ■> l b l ^ d b d ■ ■ ■ i l b l ( ^ ' ' 

This follows from contracting the Levi-Civita antisym- 
metric tensor e A k „■ u with itself and observing that 
]d°d'" l i°i _ 

the two sets of indices i and b run over d distinct values 
each. 

The Jacobian of the map from x G S d to 8 G S d C 
is 

80 



8x 



8; 



1 



— ~T, e i 6 b ■••b ®i ®b ' ' ' 8i 9 b 

d\ l d 'i °d D i L d °d *i °i 



(C39) 



f f) f) 

d\ i d - i i h d - h i l d h d 



H D i 



\d(d-l)/2 



H.b b, u ij u b 



8; 



■8, 9, 



i, u b. 



d\ l d"d'" l l"l "d "d 

The conserved quasiparticle charge then becomes 

O = [ 7° 

x£S d 



=(-i) 



<Ji .h ...i.h. 9i Oh ' ' ■ 9i 



xes d 



= H 



di 9u ■ ■ ■ <9, 



d i d b d - i i b i l d b d H b i 



x£S d 



=(-i) d+1 (-) d(d - 1)/2 rf! J N d {6) 
9es d 

(_jd(d-l)/2(_^d+l 



j e *d h d--- i i h i'^ i d b d--- i i b S °} 

0&s d 

(C40) 
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with the help of Eq. (C38) to reach the last equality. 



□ 

For convenience, we introduce another family of indices 
/i through 



Mil ' ' ' ) A*2d —1:''' 7 2d, 



Hi' ' ' ? *d — lj ' ' ' j 



(C41) 



- ■ • ■ ! U d 



rf+ 1, 



,2d. 



For any tensor J^...^ that reduces to the fully anti- 
symmetric tensor h „■ h , we have the identity 

l d°d'" l l°l 



7 =M e 



.4. 



(C42) 



since the contraction on the left-hand side of this equa- 
tion yields a combinatorial factor of (2d)! whereas the 
contraction on the right-hand side yields a combinatorial 
factor of d\ x <i!, i.e., one d\ for the family i of indices 
and another d\ for the distinct family b of indices. In 
particular, we can choose 

(C43) 

where th e sub script refers to the semi-classical Green 
function jC32T ). 

At last, we add the non-compact imaginary-time label 
with the introduction of the family v of indices, 



U ir-- -, V 2d+l =0,1,-- - > 2d : 
Mil ' ' ' i M2d = 1' ' ' ' ' 2 "> 
' ' ' 3 *d = 1 -i ' ' ' j d, 

bi, •■• ,b d = d+ 2d. 



(C44) 



Define 



:= / [ tv R (gd v g- 1 ---gd l/ g- 1 ) . 



U p£S d 



duj fdQ d (p) fdn d (e) 



K id p£S d 0£S d 



2tt J (2ir) d J (2ir) d 

s d s d 



(C45a) 



The "surface" element of S d is here denoted by dCl d . It 
follows that 



(_)^-l)/2 ( _ i)rf+ l d! 
^adia - (2d+l)j ( j 

x / e„ .... tr fl fe g- 1 ---^ a" 1 ) 

a: 

(C45b) 



thereby completing the proof of Eq. ( 3.38 ). This is the d- 
th Chern number in (2d+l)-dimensional Euclidean space 
and time. 
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